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Abstract 

Microscopic calculations show a strong parity dependence of the nuclear level density at low excitation energy of a 
nucleus. Previously, this dependence has either been neglected or only implemented in the initial and final channels 
of Hauser-Feshbach calculations. We present an indirect way to account for a full parity dependence in all steps of a 
reaction, including the one of the compound nucleus formed in a reaction. To illustrate the impact on astrophysical 
reaction rates, we present rates for neutron captures in isotopic chains of Ni and Sn. Comparing with the standard 
assumption of equipartition of both parities, we find noticeable differences in the energy regime of astrophysical 
interest caused by the parity dependence of the nuclear level density found in the compound nucleus even at sizeable 
excitation energies. 
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1. Introduction 

Nuclear reactions in systems with high level den- 
sity at low and intermediate energies are commonly 
treated in the compomid mechanism 1 2 3J. This re- 
action mechanism was first postulated by Bohr in his 
well-known independence hypothesis, stating that 
reactions can proceed via formation of a compound 
nucleus and that the decay of the compound nucleus 
is determined entirely by its energy, angular momen- 
tum, and parity, and not by the way it was formed 
[4]. This hypothesis remains valid below a projec- 
tile energy of a few tens of MeV. At higher ener- 
gies, doorway states, pre-compound, and direct reac- 
tions become increasingly important. In the follow- 
ing, we focus on the low energy region and thus on 
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the pure statistical picture (Hauser-Feshbach the- 
ory, e.g. [1I2| 1 because our ultimate goal is to propose 
an improved description for nuclear astrophysics. In 
astrophysical nuclear burning, the relevant energy in 
the projectile-target system does not exceed about 
200 keV for neutron-induced reactions and 10-12 
MeV for proton- and a-induced reactions [S]. Most 
astrophysically important reactions occur at even 
significantly lower energy. 

Angular momentum conservation is included in 
standard Hauser-Feshbach theory and thus the 
Bohr hypothesis independently holds for each spin 
J and parity tt of the compound nucleus, formed 
from the interaction of a projectile with a target 
nucleus. There are two fundamental assumptions in 
the derivation of this theory: 1) There are always 
sufficient compound-nuclear states with J^ in the 
relevant excitation energy range; 2) the wave func- 
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tions of the compound nuclear states have random 
phases, so that interferences between reactions pro- 
ceeding through different compound nuclear states 
vanish. Due to the strong energy-dependence of the 
nuclear level density, these assumptions are valid for 
most reactions (especially at intermediate energy) 
on stable targets studied in the laboratory. How- 
ever, it was shown "5] that the level density becomes 
too low for the application of the Hauser-Feshbach 
statistical model for astrophysically important re- 
actions involving nuclei far off stability, exhibiting 
small particle separation energies, or even for nuclei 
close to stability around closed shells [6] at the low 
end of astrophysically relevant energies. 

For what follows it is important that in appli- 
cations of the Hauser-Feshbach model one assumes 
that both parities are equally present in the com- 
pound nucleus at the formation energy. This pre- 
sumption clearly is not valid at very low excitation 
energy (e.g. due to pairing effects), but there are also 
many indications from theory as well as some exper- 
iments that parities may not be equilibrated even 
at considerably large excitation energies, in some 
cases up to 12 MeV |7I8I9| . Similar results are found 
in different approaches, e.g., in combinatorial meth- 
ods using single-particle energies from microscopic 
Hartree-Fock-Bogoliubov calculations [10 as well as 
in recent Shell Model Monte Carlo (SMMC) calcu- 
lations [11I12J . We note, however, that very recent 
experimental data for ^^Ni and ^"Zr are in accord 
with the equal-parity assumption at excitation en- 
ergies as low as 7 MeV (for ^^Ni in disagreement 
with the theoretical predictions). Nevertheless, the 
assumption of parity-independence of the level den- 
sity is clearly doubtful for a large number of exotic 
nuclei in the energy range important in astrophysi- 
cal environments. 

In this Letter we incorporate the parity depen- 
dence of the level density into the statistical de- 
scription of astrophysically capture reactions in all 
stages of the reaction. Parity-dependent level densi- 
ties have been used before in statistical model cal- 
culations |10|9) . But these concern the distribution 
of initial and final states, where the first can be pop- 
ulated due to the finite temperature in the astro- 
physical environment, while the parity-dependence 
in the newly formed, excited compound nucleus has 
never been considered before. Speaking in Bohr's 
terms, our modification impacts the formation cross 
section of the compound state. 

In the next section we present the details of the 
modification. This is followed by some restricted ex- 
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Fig. 1. Schematic sketcli of tlie compound capture reaction. 
A particle (neutron or proton) is captured in the state /i 
in the nucleus with mass number A, exciting a state in the 
compound/daughter nucleus {A + 1) at energy E and with 
angular momentum and parity quantum numbers J and it, 
respectively. This state decays by 7 emission to the state u 
in the same nucleus. 

amples for application to astrophysical neutron cap- 
ture which are merely given to discuss the model 
and to illustrate the possible implication for astro- 
physics. The final section gives a summary and an 
outlook to future work. 



2. Formalism 

The Hauser-Feshbach expression for the cross sec- 
tion of a reaction proceeding from the target state 
/i with spin J^* and parity n^ to a final state v with 
spin J^ and parity tt^ in the residual nucleus via a 
compound state with excitation energy E, spin J, 
and parity tt (see Fig. [1]) is given by 
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where Eij , M^ are the center-of-mass energy and 
the reduced mass, respectively, in the initial sys- 
tem, while Jj is the ground state spin of the 
projectile. The transmission coefficients T^ = 
TiE,J, tt; E^, Jt.TT^), T^ = T{E, J, tt; E^^, J^, <J 
describe the transitions from the compound state to 
the initial and final state, respectively. The sum of 
the transmission coefficients of all possible channels 
are given by Ttot- 



As indicated, the transmission coefficients depend 
on the energy E and the quantum numbers for an- 
gular momentum J and parity tt of the states ex- 
cited in the compound nucleus. They can in princi- 
ple be calculated by solving the Schrodinger equa- 
tion for the appropriate degrees of freedom. Such a 
microscopic approach will describe simultaneously 
and consistently all the states in the compound nu- 
cleus including their dependence on parity. Again 
in principle, this microscopic Schrodinger equation 
can be mapped onto a complicated optical potential 
which depends on energy and the other quantum 
numbers. There have been first attempts to derive 
such microscopic potentials, which show indeed a 
strong dependence on parity |13|14j . However, these 
attempts are not realistic enough and are restricted 
to a few scattering systems. Thus in astrophysi- 
cal (and other) applications of the Hauser-Feshbach 
model the transmission coefficients are calculated 
from optical potentials which are expected to give a 
reasonable and global account for the many nuclei 
needed in nucleosynthesis calculations. For our dis- 
cussion it is relevant that these global optical poten- 
tials do not depend on parity and hence also a pos- 
sible parity dependence of the transmission coeffi- 
cients is lost. To overcome this shortcoming we pro- 
pose here an indirect way. It is based on the observa- 
tion, that for the average transmission coefficients, 
there is the relation T oc (F) /D, involving the level 
spacing D — 1/ p and the average level width (F) in 
the considered reaction channel. This linear propor- 
tionality between transmission coefhcient and level 
density p leads us to define 
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where T{E, J) is a transmission coefficient calcu- 
lated for a global, parity-independent potential (in- 
cluding centrifugal potential) and the parity depen- 
dence is introduced by the weighting factor (tt = ±) 
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The factor 2 accounts for the proper normalization. 
This approach assumes that the parity dependence 
of the microscopic potential can be fully mapped 
onto the level density appearing in the standard 
Hauser-Feshbach equations. In the following appli- 
cations we will use the same ansatz Eq.[2]also for the 
transmission coefficients in the final channel. As all 
transmission coefficients are evaluated at the same 
energy in Eq. [H the (3 factors of the total transmis- 



sion coefficient (denominator) and one of the (3 fac- 
tors of the nominator cancel. 

As mentioned above, the Hauser-Feshbach ap- 
proach assumes a "sufficient" number of levels in 
the excited compound nucleus so that an averaged 
transmission coefficient T is justified and the model 
is applicable. For astrophysical application in the 
determination of astrophysical reaction rates the 
incident energy distribution is given by a Maxwell- 
Boltzmann distribution giving rise to a relevant 
energy window |5|3| . It has been shown that about 
10 contributing levels (depending on which partial 
waves are dominating) within this energy window 
are sufficient. This basic conclusion is not affected 
by our treatment. However, the parity dependence 
may enhance or reduce the number of available rel- 
evant levels and thus the applicability limits have 
to be reevaluated taking into account the spins and 
parities of the initial states. It should be noted that 
we do not change the total level density but just 
distribute it differently between the parities. 

3. Results 

To explore the possible effects of our modification 
we have performed a series of neutron capture cross 
section calculations. At first we have performed con- 
ventional calculations in which we assumed parity 
equipartition at all energies in target, compound nu- 
cleus and residual (calculation a). Secondly we have 
restricted the parity dependence to the level densi- 
ties of the target and residual nucleus (calculation 
b - these are similar to those of Mocelj et al. [S]). 
Thirdly, we used a parity-dependent level density in 
all three steps of the statistical treatment: the target, 
the compound nucleus and the residual nucleus (cal- 
culation c). Our calculations have been performed 
using the spin and parity dependent level densities of 
Hilaire and Goriely [10] . To explore how sensitive the 
results depend on the set of level densities adopted 
we have repeated our calculations using backshifted 
Fermi gas level densities with the parametrizations 
as derived by Rauscher and Thielemann 5j and the 
parity dependence as defined by Mocelj et al. [9] . 

In this paper we focus on (n, 7) for which the in- 
fiuence of the parity dependence of the level density 
can be discussed considering either the initial neu- 
tron capture or the final 7 decay. In the following 
we have chosen to consider the final 7 transitions 
which we assume to be either of parity-conserving 
Ml or parity-changing El multipolarity. The reac- 



tion scheme is shown in Fig.[TJ Furthermore, we show 
astrophysical reaction rates which include weighted 
sums over thermally excited states given by a ther- 
mal Maxwell-Boltzmann distribution according to 
the conditions in a stellar plasma (for the relevant 
definitions see, e.g., Ref. [3]). 

As a first example we discuss the ^*Ni(n,7)^^Ni 
reaction for which experimental data are available 
for comparison [15] . 
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Fig. 2. Stellar reaction rate of Ni(n,7) Ni; crosses: rec- 
ommended values from 'KADoNIS vO.2' 1151 . (a) without 
parity dependence, (b) using parity dependence for the final 
states, (c) using parity dependence for the final states and 
the compound formation using the level densities of |10l 

Fig.[2]shows the rate for this reaction as a function 
of temperature. One observes basically no difference 
between calculations a) and b); i.e. the consideration 
of a parity dependence in the level densities of the 
target and residual nucleus has no effect in this case. 
This is in agreement with the findings of Mocelj et 
al. [9] and is mainly caused by the fact that both cal- 
culations use the experimentally known spectrum at 
low energies. However, considering the parity depen- 
dence of the level density in the compound nucleus 
(calculation c) reduces the rate by about 30% which 
is a non-negligible effect. The origin of this reduc- 
tion becomes clear when one inspects Fig. [3] which 
shows the ratio of parity-projected level densities for 
58.59m defined as p^ = p{E, n) = J2j p{E, J, tt). We 
have summed over all spins as different values in J 
have qualitatively the same dependence in the ratio 
P-Ip+ plotted in Fig. [3l For the (n,7) reaction on 
58Ni the 7 transitions in the compound nucleus ^^Ni 
are dominated by El multipolarity at the relevant 
energies. At these energies above the neutron thresh- 
old negative-parity states dominate the spectrum of 
5^Ni due to the negative parity of the unpaired neu- 
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Fig. 3. Ratios p_/p-|_ of several nickel nuclides; the par- 
ity distribution is from Hilaire et al. |10) - we used 
p{E, n) = y^ p{E, J, 7r) here; the arrows mark the neutron 
separation energy of the corresponding nickel isotope. 

tron occupying single particle states in the pf shell. 
Thus initial states for El transitions into these states 
must have positive parity and at stellar conditions 
they have to reside just above the neutron thresh- 
old energy (S'„ = 9 MeV in ^^Ni). At such modest 
excitation energies the nuclear models predict still 
a dominance of negative parity over positive parity 
states in the level density. As we use the same total 
level density in all calculations a), b), and c), the 
ratio p-/p+ > 1 at the relevant energies yields a 
reduction of the dominant El transitions compared 
to a calculation which assumes parity equipartition, 
i.e. P-/p+ = 1. We note that this reduction gets 
smaller with increasing temperature as then higher 
excitation energies, at which the /3_/p+ gets closer 
to unity, contribute more to the stellar reaction rate. 

For Ml transitions the effect is opposite as these 
require negative-parity initial states for this reac- 
tion. These are enhanced compared to the parity- 
equipartition assumption and hence the contribu- 
tion of the Ml transitions relatively increase. This, 
however, has not much effect on the total cross sec- 
tion which is dominated by El transitions. 

It is also satisfying that our calculation yields 
a slightly better agreement with the empirical 
data from the KADoNIS compilation when parity- 
dependent level densities are incorporated into the 
statistical model (see Fig. [2]). 

Fig. [4] shows the neutron capture rates for the 
chain of nickel isotopes as obtained with full parity 
treatment (calculation c) relative to the standard 
treatment without parity dependence (calculation 
a). To understand the results one has to consider 



that the importance of parity-dependent level den- 
sities in the statistical calculation of stellar reaction 
rates depends on several ingredients: 1) the energy 
dependence of the ratio p_//3-|_, 2) the excitation en- 
ergy of the neutron threshold in the compound, 3) 
the competition of parity-changing (El) vs. parity- 
conserving (Ml) 7 transitions. 
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Fig. 4. Ratios of the stellar neutron capture rates of Ni iso- 
topes with parity influence in target, compound, and resid- 
ual to the stellar rate without any parity influence - the up- 
per graph was made by using the level densities from |10l 
and the lower graph by using the parity distribution from 
|9] combined with a level density from [5]. The ratios are 
shown for a temperature T = 1 GK of the stellar plasma. 

The ratio p_/p+ is shown in Fig. [3] for several 
nickel isotopes. Obviously positive parity states 
dominate the low-energy spectrum for even-even 
nuclei, while negative parity-states are more abun- 
dant at low energies for odd-A nuclei. With increas- 
ing energy the ratio /9_/p+ tends to unity. However, 
the energy at which parity equipartition is reached 
depends on nuclear structure, i.e. it is basically 
determined by the energy difference of the Fermi 
energy to the nearest level with different parity and 



the occupation of the levels. Thus, the equipar- 
tition is achieved at increasingly higher energies 
from ■^^Ni to ^^Ni, where the /7/2 and ^3/2 orbitals 
matter. By adding more neutrons it becomes en- 
ergetically easier to excite those to the 39/2 level 
to make levels of opposite parity. As the energy 
difference to this level decreases with increasing 
neutron number states with opposite parity (com- 
pared to the ground state parity) can be reached 
at lower energiea ^ I. However, only the excitation of 
an odd number of nucleons from or into the next 
oscillator shell changes the parity of the state. This 
leads to an oscillatory behavior in the p_/p_|_ ratios 
at energies already below the neutron thresholds 
(indicated by arrows in Fig. [3]). For ^'^Ni with the 
neutron number TV = 40 the -pj shell is completely 
occupied in the independent particle model. Hence 
parity-changing transitions appear at quite low en- 
ergies in the nickel isotopes around ^^Ni. Moving to 
even larger neutron numbers, and thus closer to po- 
tential r-process nuclei around ^*Ni, nickel isotopes 
have mainly positive parity states at low energies 
as orbitals in the gds shell have positive parity. We 
note further that equipartition is reached at some- 
what higher energies (about 3 MeV) in even-even 
nuclei than in odd-A nuclei due to pairing. 

Fig. [3] also shows the neutron separation energies 
[16j which obviously decrease with increasing neu- 
tron number along an isotope chain. The odd-even 
staggering is due to pairing. 

Usually El transitions, which are modelled by 
a Lorentzian centred around the giant dipole reso- 
nance in our approach ^, dominate over Ml tran- 
sitions, which are described by the single-particle 
model which makes the Ml strength energy inde- 
pendent [2J. However, if the capture occurs at ener- 
gies significantly below the giant dipole resonance, 
El transitions are strongly suppressed relative to 
Ml transitions in this model and the latter can dom- 
inate. This can occur in very neutron-rich isotopes 
with very low neutron thresholds. 

For the nickel isotopes with largest proton excess 
we calculate an enhancement in the stellar rates if 
the parity dependence of the level densities is con- 
sidered (Fig.|4|). For these nuclei the neutron thresh- 
olds are quite high and, due to the excitation of an 
odd number of nucleons from the ds shell, there is an 
enhancement of states in the compound with oppo- 



^ The protons do not play an important role for determining 
the parity at low energies since nickel is a closed shell nucleus 
for protons. 



site parity to the ground state. This yields a shght 
increase in the rate. The larger effect stems from the 
fact that for these nuclei low-energy states are ex- 
perimentally not known and have to be modelled. 
These states, however, are likely to have the same 
parity as the ground state (positive for even-even 
and negative for odd- A proton- rich isotopes). Con- 
sidering this fact increases the rate compared to cal- 
culations which assume parity equipartition also at 
low excitation energies. 

For the isotopes ^^"S^Ni we observe a reduction 
of the rate, if parity dependence is considered. The 
origin is the same as discussed above for the case of 
^*Ni, as El transitions are reduced as p+/ p- < 1 
{p-l p+ < 1) at the energies just above the neutron 
threshold for odd-A (even-even) isotopes. 

For the isotopes ^^-72^,^^ ^-^^ neutron thresholds 
are located at energies where the ratio of parity- 
dependend level densities is rather close to unity, 
but still shows some oscillatory behavior. As a con- 
sequence the rates are slightly enhanced or reduced 
depending on the fact whether the ratio is just above 
or below unity at the energies above the neutron 
threshold. 

For the most neutron-rich Ni isotopes, states with 
the same parity as the ground state dominate the 
spectrum at energies around the neutron threshold. 
This leads to a reduction of the rate if El tran- 
sitions dominate. However, for the capture on the 
nickel isotopes ^^Ni and ^''Ni the thresholds are so 
low (0.52 MeV and 0.17 MeV, respectively) that 
Ml transitions contribute more in our model than 
El captures; hence the rate is increased compared 
to the case where parity equipartition is assumed. 
Due to pairing the neutron threshold in ^"Ni is at 
3.7 MeV and El capture dominates. We note that 
these nuclei are close to the r-process path. Hence 
our discussion clearly shows that the effect of parity- 
dependend level densities on the neutron capture 
rate is quite sensitive to the neutron separation en- 
ergies and the competition of Ml and El transitions 
which are both not sufhciently well known yet. We 
also note that the low density of states makes the 
use of a statistical model for the very neutron-rich 
isotopes questionable. Furthermore direct neutron 
captures should contribute to the rates for these nu- 
clei. Here the parity dependence of the optical po- 
tential should be incorporated into the models and 
possible effects studied. 

We have repeated our calculation of neutron cap- 
tures on the nickel isotopes using the parity distri- 
bution of Mocelj et al. 9J combined with the back- 



shifted Fermi gas model of Rauscher et al. [5]. As 
is shown in the lower graph of Fig. [4] the effects are 
quite similar to those obtained with the Hilaire and 
Goriely level densities [10] . 

Finally we have performed a series of calculations 
for the neutron capture on the tin isotopes (Fig. [5]). 
Again El transitions dominate and for the same rea- 
sons as explained in details above (e.g. for the case 
of ^^Ni) the consideration of parity-projected level 
densities lead to a reduction of the rates. For the tin 
isotopes -'^^°~^'^°Sn the neutron intruder state ft.11/2, 
with the opposite parity to the other orbitals in the 
gds shell, plays a special role which has no equivalent 
for the nickel isotopes. It leads to a rather fast par- 
ity equilibration in the level densities which reaches 
ratios p-/ p+ close to unity at energies around the 
neutron thresholds [10 . As a consequence the rates 
for neutron captures on the mid-shell tin isotopes 
change only mildly if a parity dependence of the 
level densities is considered. For the heavier tin iso- 
topes, the TV = 82 neutron shell gap at ^'^^Sn and 
the fact that the two lowest single particle orbitals 
beyond A^ = 82 (/7/2JP3/2) have the same parity 
as the ft.11/2 intruder level have the effect that the 
equipartition of parities is reached at larger excita- 
tion energies than the respective neutron thresholds 
for tin isotopes beyond ^'^^Sn. For similar reasons as 
for the case of ^^Ni, the capture rate is reduced for 
these tin isotopes. The odd-even dependence in the 
rates are caused by the pairing effects in the neutron 
threshold energies. For the even heavier tin isotopes 
the neutron intruder state 113/2, with opposite par- 
ity to /111/2, /7/21P3/2 becomes important resulting 
in a fast parity equipartition of the level densities. 
As a consequence the capture rates do not change 
much, if parity-dependent level densities are used. 



4. Conclusion 

We have presented a simple method in the frame- 
work of the statistical Hauser-Feshbach theory to 
account for a full parity dependence including non- 
uniformly distributed parities in the nuclear level 
density of the compound nucleus. This goes beyond 
previous work which only accounted for parity- 
dependent level densities in the initial and final 
channels but not in the compound step of the reac- 
tion. We applied our method to capture reactions 
on Ni and Sn nuclei, using a parity dependence 
in all steps of the compound nucleus reaction. We 
conclude that this treatment can have a noticeable 
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Fig. 5. Ratios of the stellar neutron capture rates of Sn iso- 
topes with parity influence in target, compound, and resid- 
ual to the stellar rate without any parity influence - the up- 
per graph was made by using the level densities from 1101 
and the lower graph by using the parity distribution from 
[9] combined with a level density from [5]. The ratios are 
shown for a temperature T = 1 GK of the stellar plasma. 

effect on astrophysical reaction rates for nuclei far 
from stability. In principle, our approach can also 
be extended to include a spin dependence or a more 
general dependence on the level density of the com- 
pound nucleus. This will be the focus of future work. 
This work was performed in the framework of 
the SFB634 of the Deutsche Forschungsgemein- 
schaft and supported by the Swiss National Science 
Foundation (grant 2000-105328). 
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